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Status of jet cross sections to NNLO 
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I review the state-of-the-art for fully differential numerical NNLO programs. Topics which are covered include 
the calculation of two-loop amplitudes, multiple polylogarithms, cancellation of infra-red divergences at NNLO 
and the efficient generation of the phase space. Numerical results for e + e~ — > 2 jets are also discussed. 



o 

CO 
NO 

O 

NO 

o 

a: 

Oh! 



X 



1. INTRODUCTION 

In the near future the LHC experiments will 
provide us with a large sample of multi-particle 
final states. In order to extract information 
from this data, precise theoretical calculations are 
necessary. This implies to extend perturbative 
calculations for selected processes from next-to- 
leading order (NLO) to next-to- next-to-leading 
order (NNLO) in the perturbative expansion in 
the strong coupling constant. Due to a large vari- 
ety of interesting jet observables it is desirable not 
to perform this calculation for a specific observ- 
able, but to set up a computer program, which 
yields predictions for any infra-red safe observ- 
able relevant to the process under consideration. 

The past years have witnessed a tremendous 
progress in techniques for the computation of 
two-loop integrals and in the calculation of two- 
loop amplitudes [1]. In addition, several options 
for the cancellation of infra-red divergences have 
been discussed [2-12]. Among those, the subtrac- 
tion method - well-known from NLO computa- 
tions [13-17] - and sector decomposition [18] are 
the most promising candidates. First numerical 
results have become available for the processes 
e+e- -> 2 jets and H -> 77 [19-22]. 

2. TWO-LOOP AMPLITUDES 

The master formula to calculate an observable 
at an collider with no initial-state hadrons (e.g. 
an electron-positron collider) is given by 
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2K(s) (2Ji + 1) (2J 2 + 1) 



E 

helicity 



A n \ 2 (1) 



where p\ and p 2 are the momenta of the initial- 
state particles, 2K(s) = 2s is the flux factor 
and s = (jpx + P2) 2 is the center-of-mass energy 
squared. The factors 1/(2 Ji + 1) and 1/(2 J 2 + 1) 
correspond to an averaging over the initial helic- 
ities. dcj> n —2 is the invariant phase space measure 
for (n — 2) final state particles and Offl (pi, ...,p n ) 
is the observable, evaluated with an n-parton con- 
figuration. The index j indicates that the leading 
order contribution depends on j partons. The ob- 
servable has to be infra-red safe, in particular this 
implies that we must have in single unresolved 
limits 



o^pu-.Pn+i) - e#>(pi,...y n ), 

and in double unresolved limits 
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A n is the amplitude with n partons. At NNLO we 
need the following expansions of the amplitudes: 
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Here An denotes an amplitude with n partons 
and / loops. The master formula eq. Q is rewrit- 
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ten symbolically as 

{0 U) } = W C#> d(7 „ ( 5 ) 

n J 

and the LO, NLO and NNLO contribution as 



+ 



(o^) NNLO = J o% d*% 2 

jo%da^ 1+ JoPdaV. (6) 

To compute the two-loop amplitudes, new tech- 
niques have been developed. Examples include 
the application of the Mcllin-Barnes transforma- 
tion [23,24], differential equations [25], nested 
sums [26-29], or sector decomposition [18]. In ad- 
dition, several methods, which reduce the work- 
load have been invented. Prominent examples are 
integration-by-parts-idcntities [30,31], reduction 
algorithms [32,33] or the cut technique [34]. In 
the results, like for example in the two-loop am- 
plitude for e+e~ — > 3 jets [35,36], a new class of 
mathematical functions appears. These are the 
multiple poly logarithms. 

3. MULTIPLE POLYLOGARITHMS 

The multiple polylogarithms are defined by 
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Alternatively, they can be defined by an integral 
representation. We introduce the auxiliary func- 
tions 
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In this definition one variable is redundant due to 
the following scaling relation: 



G{z l ,...,z k ;y) = G{xz 1 , ...,xz k ;xy) 



(9) 



With the short-hand notation 

G((V^0, zi, Zk-i,0^0, z k ; y) 

mi— 1 mfc— 1 

One then finds 
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Multiple polylogarithms have interesting mathe- 
matical properties. In particular, there are two 
Hopf algebra structures defined on them [36-38] . 
At the end of the day of an analytic calculation 
physicists would like to get a number. This re- 
quires a method for the numerical evaluation of 
multiple polylogarithms. The simplest example is 
the numerical evaluation of the dilogarithm [39] : 



Li 2 (x) = -J 



dt 



ln(l - t) 
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The power series expansion can be evaluated nu- 
merically, provided |x| < 1. Using the functional 
equations 

- 2 1 
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^-i(ln(-x)) 2 , (13) 



Li 2 (x) = -Li 2 (l - x) + — ln(x) ln(l - x), 
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any argument of the dilogarithm can be mapped 
into the region |x| < 1 and —1 < Re(x) < 1/2. 
The numerical computation can be accelerated 
by using an expansion in [— ln(l — x)] and the 
Bernoulli numbers B;\ 



(-ln(l-x)) 



i+l 
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The generalisation to multiple polylogarithms 
proceeds along the same lines [40]: Using the in- 
tegral representation of 



z 2 ,...,z k ;y) 



(15) 



one transforms all arguments into a region, where 
one has a converging power series expansion. To 
accelerate the convergence of the series the Holder 
convolution is used [41]. 
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4. INFRA-RED DIVERGENCES 



The individual contributions on the r.h.s. of 
eq. © to (0^) NL0 and (0^) NNL0 are in 
general infra-red divergent, only the sum is fi- 
nite. However, these contributions live on differ- 
ent phase spaces, which prevents a naive Monte 
Carlo approach. To render the individual con- 
tributions finite, one adds and subtracts suitable 
chosen terms. The NLO contribution is given by 



+ I (o^da^+O^oda^). (16) 



The notation On ° dctn' 1 ^ is a reminder, that in 
general the approximation is a sum of terms 

0«o«fag».*> = E^^' 15 ( 17 ) 

and the mapping used to relate the n + 1 parton 
configuration to an parton configuration differs 
in general for each summand. In a similar way, 
the NNLO contribution is written as 

+ {0 (3))NNLO^ (18) 

with 

ln+2 — J (®n-\ 
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-i+2 aa n+2 w n+l ° aa n+l 



/Q(j)\NNLO 



{0^ L ° = I (O^l + O deft.? 

-OW o da^ 

+0?"i«?'"). (19) 

d&n+i i s * ne NLO subtraction term for (n + 1)- 
parton configurations, da>n' 2 ^ and dan are 
generic NNLO subtraction terms. It is convenient 
to split these terms into 

da^ = dai^L-da^ 
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such that da 
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l (o l) n are approximations to dct^+i ■ 

5. PHASE SPACE GENERATION 

It is a well-known fact, that in the collinear 
limit spin correlations remain. For example, the 
spin-dependent splitting functions for g — > qq 
reads 



p(o,i) 
g^qq 



-g 



•4z(l 
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The collinear limits occurs for k , 



— > 0. The term 



1 k^ k v 



v4/, 



(22) 



is proportional to the spin correlation. In four 
dimensions the spin-averaged splitting functions 
are obtained by integrating over the azimuthal 
angle ip of pi around p. By using spin-averaged 
antenna functions, the subtraction terms have not 
the same point- wise singular behaviour as the ma- 
trix elements, which is required for local subtrac- 
tion terms. Instead, cancellations of singulari- 
ties occurs only after an integration over the az- 
imuthal angle over all collinear splittings of the 
matrix elements. For n final-state particles, this 
is a one-dimensional integration in the (3n — 4)- 
dimensional phase space. It can be shown, that 
in the single collinear limit, the spin correlation 
depends on the azimuthal angle ip as 



1 k^ k v 



C + C 2 cos{2<p + a). (23) 



One can therefore perform the average with two 
points, where the azimuthal angle takes the val- 
ues if and ip + 7r/2, while all other coordinates 
remain fixed. In detail this is done as follows: We 
partition the phase space into different channels. 
Within one channel, the phase space is generated 
iteratively according to 



£>n+l 



d<t>nd<j)Dipole i,j,k- 



(24) 



For each channel we require that the product 
SijSjk is the smallest among all considered chan- 
nels and that sy < Sjk- Therefore it follows that 
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with channel (i,j,k) also channel (k,j,i) has to 
be included into the partioning of the phase space. 
For the dipole phase space measure we have 

1 1 2tt 

d<pdi P oie = 7^ J dy (l-j/) J dz j dip. (25) 



We can therefore generate the (n+ l)-parton con- 
figuration from the n-parton configuration by us- 
ing three random numbers u\, u 2 , uz and by set- 
ting 



2/ = tii, z = U2 if = 2ttu3. 
This defines the invariants as 

&ij ySijki ^ik ^(1 y)Sijkj 



s jk = (1 - z)(l - y)s 



ijk- 



(26) 



(27) 



From these invariants and the value of ip we 
can reconstruct the four-momenta of the [n + 1)- 
parton configuration [42]. The additional phase 
space weight due to the insertion of the (n+ l)-th 
particle is 



w = 



16tt 2 



(i-y). 



(28) 



We have therefore a parameterization of the phase 
space, such that for every collinear limit the az- 
imuthal average can be easily performed, while 
keeping all other coordinates fixed. It is clear 
that this procedure can be iterated for multiple 
collinear emissions. 

6. NUMERICAL RESULTS 

A test case for the methods discussed above 
for the cancellation of infra-red divergences and 
numerical phase space integration is the process 
e+e~ — ► 2 jets. The two-jet cross section has the 
perturbative expansion 

(a)(0) ( 1 + 2l B V-M + (^L) 2 CV 
\ 2tt \2ir / 



(29) 



(2- jet) 



(cr)(°) is the total hadronic cross section at lead- 
ing order, the numerical value at LEP energies 
v/Q 2 " = 91.187 GcV is (a)^ = 40807.4 pb. The 
jets arc defined according to the Durham jet al- 
gorithm with y = 0.01. For the numerical values 



of the NLO and NNLO coefficients B ^ et ) and 
C(2-jet) one obtains 



B {2 - Jet) = -13.674 ±0.004, 
C (2 - jet) = -231.6 ±0.3. 



(30) 



The correctness of this result can be verified 
with the help of the known results for the to- 
tal hadronic cross section at NNLO, the three-jet 
cross section at NLO and the four-jet cross sec- 
tion at LO. We have the perturbative expansions 

( CT )(*°t) = 

((T) (o) ( 1 + ^± B (tot) + (^ 1 y C (tc 

\ 2ir \2ir / 

( a )(3-jet) = 

( a )(4-^ei) = ^(0) C (4-jet)_ ^ 

The values of the coefficients are: 

B (tot) = 2j c (tot) = 5 64; ^ 

B (3-jet) = 15.679 ±0.004, 

C(3-jet) = 153 2 ± o 4) (j(4-jet) = g 4 39 ± q 

We must have 

B (tot) = B (2-jet) + _g(3-jet) j 

(j{tot) _ q(2- jet) _|_ (j(3-jet) _|_ £r(4- jet) _ ^33^ 

We find 

B (2-jet) + B (3-jet) = 2 qq 5 ± Q.006, 

C (2-jet) + C (3-jei) + ^(4-jei) = g q ± Q 5^ (34) 

which agrees well with the values of i?( tot ) and 
C (tot \ 

7. OUTLOOK 

Fully differential numerical programs at NNLO 
are a challenging task. I reviewed the state of the 
art and discussed the calculation of two-loop am- 
plitudes, the cancellation of infra-red divergences 
and methods for an efficient generation of the 
phase space. As a first application I presented 
complete results for the NNLO correction to the 
process e+e~ — > 2 jets, based on the subtraction 
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method. This gives us confidence, that the ex- 
tension to other processes like e + e~ — > 3 jets or 
pp — > 2 jets is within reach. 
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